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As the pointing accuracy and service life requirements of the DSN 70-meter 
antennas increase, it is necessary to gain a more complete understanding of the servo 
hydraulic system in order to improve system designs to meet the new requirements. 
This article develops a mat he mat i cal model for the servovalve incorporated into the 
hydraulic system of the 70-meter antenna and uses experimental data to verify the 
validity of the model and to identify the model parameters. 


I. Introduction 

JPL has long maintained the objective of improving the 
performance and operating life of the DSN 70-meter anten- 
nas. A mathematical model of the servovalve incorporated 
into the 70-meter antenna servo hydraulic system was cre- 
ated to gain a better understanding of the operation of the 
hydraulic system as one step toward this objective. During 
the development of the hydraulic servovalve model several 
approaches were considered. The first approach consisted 
of physically modeling the dynamic and steady state op- 
eration of each of the internal valve components and com- 
bining these models as a system to obtain the servovalve 
model. The internal valve components include the spool, 
bushing, feedback wire, flappers, nozzles, and torque mo- 
tor. Due to the complexity of performing this type of 
analysis and the inevitable uncertainty of the model due 
to manufacturing tolerances and complex geometries of in- 
ternal fluid flow paths, this model was determined to be 
inappropriate. An alternative model that mathematically 
describes the performance of the servovalve based upon 


experimental data was found to predict valve performance 
with reasonable accuracy. This type of model had its ad- 
vantages over the physically based model because it was 
easily implemented within standard simulation packages. 
Also, the valve manufacturer performs quality assurance 
tests to insure that the parameters of the experimentally 
based model are within allowable limits. Because of the 
stated advantages, the experimentally based model was se- 
lected to model the servovalve located within the hydraulic 
system of the 70-meter antenna. 

II. Theoretical Model Development 

A. Presentation of Mathematical Equations 

Previous experimental investigations have shown that 
the dynamic response of a flow control servovalve can be 
approximated by the second-order equations [3, 6, 7, 8] 
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The above variables, to be described in greater detail 
later in this article, are defined as 

i = input current to the valve 

i cr = critical current defining the boundary of the null 
flow region 

i H = servovalve hysteresis current 
i NB = servovalve null bias current 

is AT = lowest current input at which the flow limit of 
the valve is reached 

i T = compensating current to model threshold 
i x = X-axis intercept current adjustment 
i* x = X-axis intercept current relative to null 
K v = valve proportionality constant 
i* x — X-axis intercept current relative to null 
K v = valve proportionality constant 
I\ i = valve proportionality constant near null 
K 2 = valve proportionality constant away from null 
P c i = pressure at control port no. 1 
P c2 = pressure at control port no. 2 
Ps = supply pressure 
Pt = tank or return pressure 
Q c i — flow rate out of control port no. 1 
Q c 2 = flow rate out of control port no. 2 
t = time 

U = effective command signal 
p ~ fluid density 
w n = valve natural frequency 
£ = damping ratio 

B. Pressure-Drop-to-Flow Relationship 

A servovalve operates by metering fluid flow from the 
supply line to one control port and from the other con- 
trol port to the return line by using a restriction with a 
pressure-drop-to-flow relationship equivalent to that of an 
orifice. Using the orifice model, the fluid flow rate from the 
supply line to the control ports and from the control ports 
to the return line varies with the square root of the pres- 
sure difference across the valve and as the square root of 
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the inverse of the fluid density. Equations (3) through (6) 
provide a mathematical description of the orifice behavior. 

C. Flow Gain 

Servovalves are manufactured to produce a linear flow 
gain with respect to input current for a given valve pres- 
sure drop. In practice the flow gain is modeled as a piece- 
wise linear function. There are three linear regions of the 
flow curve as shown in Fig. 1. The three regions are re- 
ferred to as the null region, the normal flow region, and the 
flow saturation region. Within the null region where the 
current input to the valve is ±i cr of the null position, the 
valve exhibits a linear flow gain K\ according to Eq. (9). 
The value of the flow gain I\\ is dependent upon the lap 
of the valve. Critically lapped valves exhibit a flow gain 
that is equivalent to the flow gain AT 2 found in the region 
of normal flow while valves that are overlapped and un- 
derlapped have flow gains that are lower and higher than 
I\ 2 , respectively. The region of normal flow is the designed 
operating range of the valve where the rated flow gain K 2 
occurs. Once the current input to the valve exceeds a cer- 
tain level, flow saturation occurs. When saturation exists 
the flow gain decreases in a nonlinear fashion until the flow 
limit of the valve is reached. At the flow limit of the valve 
an increase in current input to the valve results in no in- 
crease of flow to the control port. The behavior of the 
flow gain in the flow saturation region is governed by such 
things as the stability of the supply pressure at high valve 
flow rates and the mechanical limitations placed on the 
valve spool. Since the supply pressure to the servovalve 
is fairly constant on the 70-meter antenna, the region of 
saturation is modeled as an abrupt change from the region 
of normal flow to the valve flow limit. 

D. Null Bias 

The null bias of a valve is the input current required to 
bring the valve to its null position under a certain set of 
operating conditions not including the effects of hysteresis. 
The null bias of the valve is modeled in Eq. (7) by the 
variable ijvB* 

E. Threshold 

Threshold is the increment of current input required to 
produce a change in valve flow rate when changing the di- 
rection of the applied current to the valve. A valve flow 
curve illustrating threshold is shown in Fig. 2. Threshold 
is mathematically modeled in Eq. (7) by the variable i? 
and is computed once every sampling period that a cur- 
rent signal is sent to the servovalve by the following logic 
statements: 
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else if ik < iB then 
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(33) 

where 


Fk = logical variable that indicates the direction of 
application of the input current. The variable 
takes on the values of —1 for decreasing current 
and *f 1 for increasing current. 

i = input current to the valve 

= upper endpoint of current range where flow does 
not change 

i B = lower endpoint of current range where flow does 
not change 

ik = input current to the valve at sample interval k 
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ix = input current adjustment required to model 
threshold 

= threshold current 
Jc = sample interval 

e = amplitude of input current signal noise 

Equations (13) through (20) are designed to keep track 
of whether the input current is increasing or decreasing 
while taking into consideration the amplitude of the input 
signal noise. A change in the direction of current applica- 
tion is detected by Eq. (21), while Eqs. (22) through (27) 
are used to redefine the current region where invariant flow 
occurs, as shown in Fig. 3. The adjustment to the input 
current required to model the threshold phenomenon is 
computed using Eqs. (28) through (32). Equation (33) is 
used to increment the sampling interval index. The logical 
statements used to model threshold have been restated in 
flow diagram format in Fig. 4. 

F. Hysteresis 

During a single cycle of the input signal, a given flow 
rate is achieved twice: once as the signal increases and once 
as it decreases. Hysteresis is the difference in valve current 
inputs required to produce this flow rate. The magnitude 
of hysteresis increases with the amplitude at which the 
input signal is cycled, as shown in Fig. 5. When the cur- 
rent is cycled at low amplitudes, the hysteresis effects are 
contained within the threshold model. As the amplitude 
of the current cycle is increased, additional modeling is 
required to predict valve performance. During the opera- 
tion of the hydraulic system of the 70-meter antenna, the 
servovalve will operate closely around a single input cur- 
rent as opposed to cycling through its full operating range. 
Under these conditions the threshold model is sufficient to 
describe the hysteresis effects of the valve. Because the 
slewing of the antenna operates over the full region of the 
valve, the flow curve of the valve may shift within the 
limits of the valve hysteresis depending upon the current 
signal history prior to antenna tracking. The shift of the 
flow gain curve due to hysteresis is modeled by the variable 
in in Eq. (7). The variable iff is treated as a statistical 
variable that changes over extended periods of time and 
takes on values within a specific distribution. 

III. Experimental Investigation 

A. Overview 

During the investigation of the performance character- 
istics of the Moog 72-163B servovalve at JPL, experimen- 
tal programs were undertaken both at Moog under the 


direction of Peter Hames and at Fluid Technologies, Inc., 
under the direction of Elizabeth Carrell. The results gath- 
ered during these investigations were used to identify the 
parameters of the theoretical model presented in the pre- 
vious section and to estimate the accuracy of the model. 

B. Frequency Response 

The frequency response characteristics of the Moog 
72-163B servovalve were determined through experiments 
conducted at Moog on January 27, 1988. It was concluded 
from an analysis of the results that the dynamic response 
of the servovalve is best approximated by a second-order 
model with a natural frequency of 55 hertz and a damp- 
ing ratio of 0.8. Both the experimental frequency response 
and the second-order model approximations are presented 
in Figs. 6 and 7. 

C. Hysteresis and Threshold 

The hysteresis and threshold effects of the Moog 72- 
163B servovalve no. 110 were experimentally determined 
at Moog on April 19, 1988. The experimental results in- 
dicated a hysteresis of 0.1 milliampere when the valve was 
cycled between negative and positive full flow. A thresh- 
old current of 0.01 milliampere was also observed during 
the testing of the valve. 

D. Steady State Flow Characteristics 

An experimental program was conducted at Fluid 
Technologies, Inc., to determine the effect that the geomet- 
rical dimensions of the valve spool have upon servovalve 
flow characteristics. The program was conducted between 
December 31, 1990, and January 7, 1991. The spool di- 
mensions along with the supply pressures for each of the 
experiments are shown in Fig. 8. Details of the experimen- 
tal hardware and procedures are presented in [1,2, 4, 5]. 

E. Data Reduction 

The valve flow constant K v and the X-axis intercept 
of the flow curve were computed by performing a least 
squares fit of the experimental flow rate versus input cur- 
rent data using the computer program MATLAB. The 
data points contained in the negative, null, and positive 
flow regions of the servovalve were fit to a first-order poly- 
nomial separately in order to conform to the piecewise lin- 
ear theoretical model. Some of the data points obtained at 
the higher current input magnitudes were removed prior 
to fitting the model to the data. The removal of the data 
points believed to be invalid was done because the valve 
was observed to saturate at flow rates near the maximum 
output capability of the pump supplying fluid to the ser- 
vovalve. Therefore, it was suspected that the supply pres- 
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output capability of the pump supplying fluid to the ser- 
vovalve. Therefore, it was suspected that the supply pres- 
sure to the valve was not constant at the larger input cur- 
rent magnitudes. There was no way of knowing whether 
the supply pressure was actually constant since the supply 
pressure was not monitored continuously during the ex- 
periments. Hence, some data points were removed where 
saturation of the valve was observed. The parameters that 
were identified from the least squares fit of the data are 
presented in Table 1, while a plot of both the experimental 
and theoretical valve flow curves are shown by the exam- 
ple in Fig. 9. It is evident from an examination of Fig. 9 
that the experimental flow curves of the positive and neg- 
ative flow regions did not necessarily match the flow curve 
of the null region from a continuity standpoint. This lack 
of continuity occurs because the data taken in the null re- 
gion were obtained in a separate experimental run from the 
data taken in the normal flow region. Hence, the observed 
discontinuity is attributed to the effects of hysteresis and 
threshold. For this reason, the data taken in the null re- 
gion were used only to identify the valve flow constant. 
The results presented in Table 1 indicate that the valve- 
flow-gain constant, K v , is slightly different for the positive 
and negative flow regions contrary to expectation. This 
small difference is believed to be the result of the valve 
spool and bushing not being “exactly” symmetrical about 
null. When modeling the valve in computer simulation 
programs, an average of the valve flow constant in the posi- 
tive and negative flow regions should be used to model the 
valve in the region of normal flow. The null bias of the 
valve was taken to be the midpoint of where the positive 
and negative flow curves intersect the X-axis. The value 


of the X-axis intercept current, used in the servovalve 
model is taken to be one half the current magnitude be- 
tween the locations where the positive and negative flow 
curves intersect the X-axis. The current value for the vari- 
able i cr that defines the boundary of the null flow region is 
computed by finding the current at which the linear flow 
curve of the null region intersects the linear flow curve of 
the normal flow region. The servovalve model parameters 
for each of the experimental valves is presented in Table 2. 

F. Summary of Results 

Experimental methods were implemented to determine 
the parameters for the Moog 72-163B servovalve according 
to the servovalve model presented in Section II. The model 
parameter values that should be used when simulating the 
valve in the hydraulic system of the 70-meter antenna are 
presented in Tables 2 and 3. 


IV. Conclusion 

In order to more completely understand the dynamic 
behavior of the servo hydraulic system of the 70- meter an- 
tennas located in the DSN, a mathematical model for a 
servovalve was developed. Experimental data were used 
to identify the model parameters for the Moog 72-163B 
servovalves that are incorporated into the hydraulic sys- 
tem of the 70-meter antenna. The servovalve performance 
predicted by the theoretical servovalve model was found 
to be in excellent agreement with the experimental data. 
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Table 1. Servovalve parameters. 


Negative flow region Null flow region Positive flow region 


Valve ID 

K V1 

m 2 /mA 

Ix, 

mA 

(m 3 /sec) 2 

K vt 

m 2 /mA 

Ix, 

mA 

^ 2 , 

(m 3 /sec) 2 

K VJ 

m 2 /mA 

Ix, 

mA 

^ 2 » 

(m 3 /sec) 2 

A-l 

1.08E-06 

— 5.43E— 01 

1.89E-10 

3.90E-07 

— 3.94E— 01 

8.30E-10 

1.18E-06 

3.01E— 01 

1.40E-10 

A-2 

1.05E-06 

— 3.73E— 01 

2.45E-10 

2.83E-06 

2.43E-01 

6.81E-10 

1.13E-06 

6.77E-01 

6.86E-11 

A-3 

1.08E-06 

— 1.46E-01 

2.53E-10 

1.99E-06 

5.85E-01 

5.63E-10 

1.07E-06 

8.87E-01 

1.03E-09 

A-4 

9.99E-07 

— 6.17E— 01 

5.40E-11 

1.03E-06 

8.86E-02 

1.71E-09 

1.04E-06 

5.29E-01 

8.30E-11 

B-l 

1.15E-06 

— 1.42E+00 

4.51E-11 

7.90E-08 

— 8.68E— 01 

5.60E-12 

1.20E-06 

2.80E-01 

1.96E-10 

B-2 

1.06E-06 

— 4.53E-01 

7.88E-11 

3.38E-07 

-1.32E-01 

2.73E-11 

1.11E-06 

5.74E-01 

1.06E-10 

B-3 

1.17E-06 

— 4.29E— 01 

7.52E-11 

1.45E-07 

3.28E-02 

3.12E-12 

1.23E-06 

4.36E-01 

5.55E-11 

B-4 

1.02E-06 

— 3.16E-01 

8.06E-11 

8.72E-07 

4.21E-02 

2.S2E-10 

1.05E-06 

2.18E-01 

9.37E-11 

B-5 

1.13E-06 

— 3.22E-01 

1.57E-10 

1.01E-06 

— 3.63E— 03 

1.80E-10 

1.17E-06 

1.15E-01 

1.43E-10 

C-l 

1.05E— 06 

-6.60E— 01 

5.35E-11 

3.27E-07 

— 2.10E— 01 

9.29E-11 

1.15E-06 

4.32E-01 

7.28E-11 

C-2 

1.17E-06 

— 8.48E-01 

1.60E-10 

5.11E-07 

— 6.99E— 01 

2.27E-11 

1.26E-06 

— 2.07E— 01 

6.94E-11 

C-3 

1.04E-06 

-7.54E-01 

8.08E-11 

1.13E-07 

— 8.64E— 02 

5.00E-12 

1.10E-06 

6.58E-01 

1.34E-10 

C-4 

1.02E-06 

— 7.96E— 01 

4.36E-11 

3.00E-09 

3.24E-01 

5.00E-13 

1.06E-06 

1.15E+00 

1.46E-10 

D-l 

1.46E-06 

7.02E-03 

4.41E-10 

1.08E-06 

4.12E-02 

4.91E-11 

1.53E-06 

3.10E— 01 

2.76E-10 

D-2 

1.23E-06 

— 5.14E— 01 

1.95E-10 

1.14E-06 

3.02E-01 

1.32E-10 

1.38E-06 

6.73E-01 

9.01E-10 


K v — valve proportionality constant. 
I x = X-axis intercept. 
a 2 = variance of the residuals. 
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Table 2. Servovalve model parameters. 


Valve ID 

Ki, 

m 2 /mA 

k 2 , 

m 2 /mA 

‘x. 

mA 

*cr» 

mA 

mA 

A-l 

3.89E-07 

1.13E-07 

4.22E— 01 

6.45E-01 

— 1.22E— 01 

A- 2 

2.83E-06 

1.09E-06 

5.25E-01 

— 3.27E— 01 

1.52E-01 

A-3 

1.99E-06 

1.07E-06 

5.17E-01 

— 6.08E— 01 

3.71E-01 

A-4 

1.03E-06 

1.02E-06 

5.73E-01 

— 6.85E+01 

— 4.39E— 02 

B-l 

7.91E-08 

1.18E-06 

8.52E-01 

9.13E-01 

— 5.71E— 01 

B-2 

3.38E-07 

1.09E-06 

5.14E-01 

7.45E-01 

6.04E-02 

B-3 

1.45E-07 

1.20E-06 

4.32E-01 

4.92E— 01 

3.08E-03 

B-4 

8.72E-07 

1.04E-06 

2.67E-01 

1.69E+00 

— 4.89E— 02 

B-5 

1.01E-06 

1.15E-06 

2.19E-01 

1.80E+00 

— 1.04E— 01 

C-l 

3.27E-07 

1.10E-06 

5.46E-01 

7.76E-01 

-1.14E-01 

C-2 

5.11E-07 

1.22E-06 

3.21E— 01 

5.53E-01 

— 5.27E— 01 

C-3 

1.13E-07 

1.07E-06 

7.06E-01 

7.89E-01 

— 4.80E— 02 

C-4 

2.85E-09 

1.04E-06 

9.72E-01 

9.75E-01 

1.77E-01 

D-l 

1.08E-06 

1.49E-06 

1.51E-01 

5.48E-01 

1.58E-01 

D-2 

1.14E-06 

1.31E-06 

5.93E-01 

4.65E+00 

7.96E-02 


K\ = valve proportionality constant in the null flow region. 

K 2 = valve proportionality constant in the normal flow region. 

i* = X-axis intercept relative to null. 

t cr = critical current defining the null flow region. 

iN B — null bias. 


Table 3. Servovalve model parameters. 


Parameter 

Value 

ijj, mA 

±0.1 

is at, mA 

20 (rated current) 

mA 

0.01 

Ps, bar 

172.4 

Pt , psi 

0 

p f kg/m 3 

854 

Wni Hz 

55 

£ 

0.8 


229 



CURRENT 



Fig. 1. Servovalve flow curve at a constant pressure drop. 


*C7 

Fig. 2. Threshold loops. 


FLOW 



Fig. 3. Invariant flow region. 



FLOW 



Fig. 4. Threshold model. 
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- EXPERIMENTAL RESULTS 
• SECOND ORDER APPROXIMATION 


EXPERIMENTAL RESULTS 

SECOND ORDER APPROXIMATION 


| ^ 

CL 


FREQUENCY, Hz 


FREQUENCY, Hz 


Fig. 6. Magnitude frequency response of Moog 72-163B servo- Fig. 7. Phase angle frequency response of Moog 72-163B servo- 


valve; supply pressure = 172.4 bar. 


valve; supply pressure s 172.4 bar. 
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VALVE 

ID 

SPOOL DIMENSIONS 

PRESSURE 

(PSI) 

A 

(In.) 

B 

(in.) 


D 

On.) 

E 

(mm) 



H 

(mm) 

A-1 

0.999 

1.003 

1.001 

0.74975 

0.03 

wm 

0.03 

0.03 

2500 

msm 

?f 

” 



0.08 


0.09 

0.09 

2500 

a 

tl 

” 

II 


0.15 


0.16 

0.15 

2500 

MSM 


tl 


99 

0.32 


0.34 

0.33 

2500 



2500 

2500 

1800 

2500 

1800 

2500 

1800 

2500 

2500 

1800 



0.9975 0.997 0.998 0.7499 0.05 0.05 0.06 0.14 


0.005 AXIAL SCRATCH EQUALLY SPACED THREE PLACES 
0.010 AXIAL SCRATCH EQUALLY SPACED THREE PLACES 
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